Orbital selective phase transition induced by different magnetic states: 
A dynamical cluster approximation study 
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Motivated by the unexplored complexity of phases present in the multiorbital Hubbard model, 
we analyze in this work the behavior of a degenerate two-orbital anisotropic Hubbard model at half 
filling where both orbitals have equal bandwidths and one orbital is constrained to be paramagnetic 
(PM), while the second one is allowed to have an antiferromagnetic (AF) solution. Such a model 
may be relevant for a large class of correlated materials with competing magnetic states in different 
orbitals like the recently discovered Fe-based superconductors. Using a dynamical cluster approxi- 
mation we observe that novel orbital selective phase transitions appear regardless of the strength of 
the Ising Hund's rule coupling J z . Moreover, the PM orbital undergoes a transition from a Fermi 
liquid (FL) to a Mott insulator through a non-FL phase while the AF orbital shows a transition 
from a FL to an AF insulator through an AF metallic phase. We discuss the implications of the 
results in the context of the Fe-based superconductors. 
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PACS numbers: 71.10.Fd,71.27.+a,71.30.+h,71.10.Hf 

The existence of an orbital selective phase transition 
(OSPT) induced by the interplay of a narrow band of 
localized electrons and a wide band of itinerant elec- 
trons has evoked considerable interest since the discov- 
ery of exotic phase transitions in Ca2- :c Sr, c Ru.04 [HQ- 
Presently, OSPT are being intensively investigated as an 
alternative to conventional Mott transitions happening 
simultaneously in all orbitals in correlated systems [Si- 
ll] . Previous works on OSPT were either based on de- 
generate orbitals with inequivalent bandwidths at inte- 
ger filling [3-9] or focused on orbitals with degeneracies 
lifted completely or partially by the crystal field split- 
ting at different fillings [13, El- Recently a new class 
of multiorbital systems -the Fe-based superconductors- 
where superconductivity appears to be strongly related 
to magnetism has opened a new field of debate. In the 
present work, we investigate a minimal two-band Hub- 
bard model of relevance to understand the nature of the 
antiferromagnetic metallic state and the experimentally 
observed reduced magnetic moment in Fe-based super- 
conductors. We will show that inclusion of (short-range) 
spatial fluctuations which are present in the dynamical 
cluster approximation (DCA) allows also for a possible 
non-Fermi liquid phase. Such a phase has been recently 
discussed for BaFe 2 (Asi_ a; P a: )2 [H and FeSe [j~3]. Our 
results suggest also a new mechanism for an OSPT in- 
duced by different magnetic states in different orbitals. 

We consider a degenerate two-orbital anisotropic Hub- 
bard model at half filling where both orbitals have equal 
bandwidths and one orbital is constrained to be paramag- 
netic (PM) (PM orbital), while the second one is allowed 
to have an antiferromagnetic (AF) solution (AF orbital). 
This model should capture the coexistence of frustrated 
bands where magnetic order is suppressed with unfrus- 




FIG. 1. (Color online) The U — T phase diagram of (a) AF 
and (b) PM orbitals for J z = [7/8 and U' = 3(7/4. The non- 
FL and AF metal are determined by the self-energy neither 
showing divergent behavior nor approaching zero as the Mat- 
subara frequency goes to zero and the spin-dependent sym- 
metry breaking field remaining finite in the metallic state, 
respectively. 



trated bands which can order antifcrromagnetically due 
to perfect nesting, as observed in bandstructure calcula- 
tions 14 1 . We employ the DCA [l| with a cluster size 
of N c = 4 in combination with a weak-coupling contin- 
uous time quantum Monte Carlo algorithm [16| . Such 
an approach ensures the consideration of the dynamical 
fluctuations and the inclusion of a symmetry-breaking 
state as well as spatial fluctuations which cannot be cap- 
tured in most of previous OSPT works usin g a s in g le-site 



dynamical mean field theory (DMFT) 

We find a novel OSPT: In the PM orbital, a phase 
transition from a Fermi liquid (FL) to a non-FL state 
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followed by a Mott metal-to-insulator transition (MIT) 
is detected, while in the AF orbital two phase transitions 
from a PM metal to an AF insulator through an AF 
metallic state are observed, regardless of the strength of 
the Ising-type Hund's rule coupling J z . All the transi- 
tions occur at different critical values of U/t. Note that 
here the OSPT occurs when the PM orbital undergoes a 
Mott MIT while the AF orbital remains an AF insulator. 
Furthermore, we also discuss the nature of the Mott MIT 
in the PM orbital at J z = where the ferromagnetic in- 
teraction between PM and AF orbitals vanishes. Such a 
consideration simultaneously allows for the study of (i) 
strong orbital fluctuations between PM and AF orbitals 
as well as (ii) the effect that the AF order in the AF or- 
bital induces on the PM orbital. Moreover, comparable 
to the single-orbital Hubbard model within a small clus- 
the cooperation between the spatial AF order 
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ter 

due to perfect nesting and frustration induced by short- 
range spatial fluctuations remains in our model. Never- 
theless, while the single-orbital Hubbard model shows a 
first-order Mott MIT 17,ll8|, we observe a continuous 



Mott MIT in the PM orbital with a critical interaction 
U c /t = 4.0 close to that obtained for the single-orbital 
system by the same DCA method with a cluster size of 
N c = 4 18, 2(|. Finally, we discuss the implications of 



FL, non-FL and Mott insulator phases are present in the 
PM orbital. The metallic regions are shrunk due to the 
enhancement of spatial AF correlations with decreasing 
temperature. The main difference between single-orbital 
and our two-orbital Hubbard model is that the first-order 
MIT present in the single-orbital Hubbard model, is re- 
placed by a continuous MIT in the two-orbital Hubbard 
model (see Fig. Q] (b)). The coupling between the PM 
orbital with spatial AF correlations and the AF orbital 
with AF order suppresses the first-order transition. Fur- 
thermore, in the weak-coupling region, while the AF in- 
sulator is present in the cluster-DMFT calculation of the 
single-orbital Hubbard model [22| , in our AF orbital the 
AF symmetry breaking field is completely suppressed due 
to thermal and orbital fluctuations between PM and AF 
orbitals (see Fig. Q](a)). In fact, in the intermediate re- 
gion the AF orbital shows an AF metallic phase and as 
the interaction U /t is increased, the AF insulator state 
is reached. The orbital selective phase where a metallic 
state in the PM orbital and an insulating state in the AF 
orbital coexist, is clearly observed in the intermediate 
regime. Such a phase is induced by different magnetic 
states in the two orbitals. This mechanism is distinct 
from previously proposed ones, such as change of filling in 



the non-Fermi liquid phase as well as the small magnetic 
moment in the AF orbital in the context of the Fe-based 
superconductors. 

The degenerate two-orbital anisotropic Hubbard model 
with a PM and an AF orbital at half-filling on the square 
lattice can be written as 

(ij)ma im 

+ J~]( U> ~ S a-a'Jz)nil a ni2a', (1) 
iaa' 

where Ci m(T (c| mo .) is the annihilation (creation) operator 
of an electron with spin a at the i-th site with orbital 
index m — (1,2). t m is the hopping matrix element 
between site i and j, U and U' are intra-orbital and 
inter-orbital Coulomb repulsion integrals, respectively, 
and J z ni\ a ni2a is the Ising-type Hund's rule coupling 
term [Hj]. The bandwidth is W = 8 (t = 1). In order to 
simulate our anisotropic Hubbard model with a PM and 
an AF orbital, a spin-dependent symmetry breaking field 
is applied on the AF orbital in the first iteration, while 
we always keep the PM solution under the condition of 
GpM{iu n ) = \{G^{iu n ) + Gi(iuj n ) for the PM orbital 
at each iteration. The converged Green's functions are 
obtained from the DCA and Dyson equations. 

First, we would like to discuss the phase diagram for 
J z = U/8 and U' = 3C//4 in the U - T plane shown 
in Fig. [TJ Analogous to the phase diagram of the single- 
orbital Hubbard model with a PM solution obtained from 



B, EI, 



non-degenerate or partially degenerate orbitals 
variation of bandwidth in different orbitals or dif- 

ferences of frustration strength in different orbitals [§]. 

In the following, we will analyze the various states in 
our phase diagram. Figs. [2] (a)-(d) show the density of 
states (DOS) for the two orbitals for several values of the 
interaction strengths U/t with U' = U/2 and J z = U /A 
at T/t = 1/24, obtained using the maximum entropy 
method. In the weak-coupling region (U/t = 0.8) both 
the PM and the AF orbital are in the PM metallic phase. 
This results in the same DOS for both orbitals includ- 
ing a quasiparticle peak at the Fermi level as shown in 
Fig. [2] (a). As the interaction increases (U/t = 1.2), the 
AF orbital - while remaining metallic - starts to display 
AF order as shown in Fig. [5] (b) . In this case, the spin up 
DOS is pairwise equal to the spin down DOS on each site. 
The PM orbital remains in the PM metallic phase with a 
quasiparticle peak at the Fermi level. In the intermediate 
region at U/t = 2.6 (Fig. [2] (c)) the metallic PM orbital 
coexists with the insulating AF orbital. When the inter- 
action is further increased to the strong-coupling region 
(U/t = 3.4), insulating states appear in both orbitals as 
displayed in Fig. [2] (d). Note that we found the OSPT 
behavior for all investigated values of J z . 

Next, we would like to clarify the nature of the Mott 
MIT in the PM orbital. In the case of the single-orbital 
Hubbard model, the double occupancy shows hystere- 
sis with coexistence regions indicating a first-order phase 
transition [l7j . Since our two-orbital system allows for 
additional orbital degrees of freedom, it is interesting to 
check whether the nature of the transition is changed due 
cluster-DMFT with N c — 4 [17] ■ we observe that the to the coupling between AF and PM orbitals. In Fig. [31 
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FIG. 2. (Color online) Density of states p(lu) per site for 
different values of the interaction strength U /t, with U' = U/2 
and J z = U/4 at T/t = 1/24. Here, only the p(ui) for spin up 
on one sublattice is shown when the AF orbital is in an AF 
state. 




FIG. 3. (Color online) The double occupancy of the PM or- 
bital as a function of U/t for three values of U' and J z at 
T/t — 1/24, with error bars as indicated. Inset: Quasiparti- 
cle weight of the PM orbital s a function of U/t. 



we plot the double occupancy as a function of U /t with 
different V and J z at T/t = 1/24. While U c /t increases 
from 3.0 to 3.5 and 3.6 when J z is decreased from U/4 
to U/8 and to with the constraint of U = U' + 2J Z , in 
none of the cases a signature of hysteresis in the double 
occupancy is found. The disappearance of the first-order 
transition is attributed to (i) the interaction between the 
orbitals and (ii) the cooperation of spatial AF correla- 
tions in the PM orbital and AF order in the AF orbital 
at all values of J z . In the inset of Fig. [3] we also show the 
evolution of the spectral weight for the PM orbital. 
Recently, the non-FL behavior in the multi-orbital as 
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FIG. 4. (Color online) Self-energies at the Fermi surface K = 
(tv, 0) as a function of (a) Matsubara frequency and of (b)-(d) 
real frequency oj/t at (b) U/t = 0.8, (c) U/t = 2.6 and (d) 
U/t = 3.4 for U' = U/2, J z = U/4 at T/t = 1/24 on the PM 
orbital. 



well as in the single-orbital Hubbard model has been 
extensivel y st udied within DMFT or its cluster exten- 
sion ^SlilllMSIiS-il. Also, non-FL behavior has 
been observed in the multiorbital Fe-based superconduc- 
tors [12]. Using the maximum entropy method [2rj| for 
the analytical continuation of the self-energy, we analyze 
the self-energy £(a>) at the Fermi surface K = (it, 0) of 
the PM orbital as a function of the real frequency oj (see 
Figs. [1(b)- (d)) for V = U/2, J z = U/4 at T/t = 1/24. 
Im Y,(iuj n ) as a function of Matsubara frequency is also 
shown in Fig. |4](a) In Fig.|4](b), ImS(cj = 0) approaches 
zero indicating a FL behavior. Upon increasing the in- 
teraction in Fig. [4] (c), the system still shows metallic 
behavior but Im£(cj = 0) has a finite value (see also 
Fig. [4] (a)). This suggests a non-FL behavior In 
the strong-coupling region (U/t = 3.4), the self-energy 
ImH(w — 0) diverges (Fig. H (d)) and the system is in 
a Mott insulating state. Also shown in Figs. 2] (b)-(d) 
is the ReS(w) obtained from the Kramers-Kronig rela- 
tion. From these results we infer that the gap opening 
in the PM orbital is caused by a blocking of the electron 
delocalization due to the Coulomb repulsion (Mott insu- 
lator), while in the AF orbital the quasi-particle peak is 
split into two parts above and below the Fermi level due 
to the antiferromagnetism with a gap of A w mU (AF 
insulator), where m is the magnetization and U is the 
interaction strength. Here we would like to stress that 
the gap opening in the AF orbital - which favors antifer- 
romagnetism - happens at a lower interaction strength 
than it happens in the PM orbital. Therefore, a Mott 
MIT occurs in the PM orbital while the AF orbital re- 
mains an AF insulator. As mentioned above, this indi- 
cates a novel OSPT driven by different magnetic states 
in different orbitals. 
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FIG. 5. (Color online) The staggered magnetization of the 
AF orbital as a function of U/t for three values of U' and Jz. 



Let us now analyze the AF metallic phase in the 
AF orbital at intermediate couplings. Unlike the high- 
T c cuprates with an AF insulating state in the normal 
phase, the undoped Fe-based superconductors show AF 
metallic behavior with a small ordered magnetic moment. 
DFT calculations [28|-|30j overestimate the values of the 
ordered magnetic moment compared to those obtained 
from experiments. With the consideration of an arti- 
ficially negative U in the DFT calculations, the small 



magnetic moment is recovered [31|, |32|. Nevertheless 



the mechanism for the reduced ordered magnetic mo- 
ment still remains controversial [33l - l38j . Very recently, 
the authors proposed a mechanism of coupling between 
frustrated and unfrustrated orbitals within a two-orbital 
Hubbard model, solved by sing le-site DMFT [§] as a pos- 
sible mechanism for the reduced magnetic moment and 
the AF metallic state. In the present work we study 
a related model (Eq. (HJ)) by DC A where spatial fluc- 
tuations - absent in the single-site DMFT calculation 
- are partially considered. In the present model, keep- 
ing one orbital always in the PM state can be viewed as 
an effective frustration due to the suppression of long- 
range correlations by the small cluster size 17j ■ In Fig. [5] 
we present the staggered magnetization in the AF or- 
bital as a function of U/t with different V and J z for 
T/t = 1/24. We observe that, as the interaction strength 
is increased, the staggered magnetization increases con- 
tinuously for all values of J z and remains small in the 
range U/t — 1.5 — 2.5 relevant for the Fe-based supercon- 
ductors. This is also consistent with DMFT results @. 

In conclusion, we have studied the two-orbital Hubbard 
model at half filling where one orbital is constrained to be 
in the PM state and the second orbital is allowed to have 
an AF solution. Analysis of double occupancy, magneti- 
zation and DOS shows that as U /t increases, while a con- 
tinuous Mott MIT through non-FL behavior is observed 



in the PM orbital, two phase transitions (PM metal to 
AF metal and AF metal to AF insulator) are detected 
in the AF orbital. We find a novel OSPT where a MIT 
occurs in the PM orbital while the AF orbital remains 
an AF insulator. We ascribe the OSPT detected in this 
work to the different magnetic states in the different or- 
bitals since other effects such as different bandwidths, 
crystal field splitting and change of band filling have been 
avoided. Even with inclusion of spatial fluctuations, the 
AF metal still survives in the AF orbital in a certain 
range of interaction strengths. These results support 
the scenario of a reduced ordered magnetic moment and 
the existence of an AF metallic state in the Fe pnictides 
driven by coupling between frustrated and unfrustrated 
orbitals [9j . Finally, by investigating the self-energy as a 
function of Matsubara and real frequency u/t, we iden- 
tify an interaction regime where non-FL behavior can be 
observed in the PM orbital. This could be of relevance for 
BaFe2(Asi_ x P a: )2 11211 and FeSe where non-FL behavior 
has been discussed jl 31 ] - 
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